Search for n° — > z/ M z/ M Decay in LSND 
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We observe a net beam-excess of 8.7±6.3 (stat) ±2.4 (syst) events, above 160 MeV, resulting from 
the charged-current reaction of and/or P M on C and H in the LSND detector. No beam related 
muon background is expected in this energy regime. Within an analysis framework of it — > v^v^, 
we set a direct upper limit for this branching ratio of r(7r° -^v^v^) /r(7r° — >aZZ) < 1.6 x 10~ 6 at 90% 
confidence level. 



The observation of the decay 7r° — > vv would imply 
new interesting physics. The pion has zero spin and 
odd intrinsic parity (i.e. J p = 0~), and it is repre- 
sented by a wave function which has the space trans- 
formation properties under inversion and rotation of a 
pseudoscalar. Momentum and angular momentum con- 
servation require that the decay v and v possess the 
same helicity. This decay provides an ideal laboratory 
to search for the pseudoscalar (P) weak interaction, be- 
cause only the P interaction allows massless neutrinos 
and antineutrinos with the same helicity in the final 
state. Furthermore, if the neutrino mass is not zero and 
the Z a couples to the right-handed neutrino with stan- 
dard weak-interaction strength, the branching ratio (BR) 
B(ir° — > v^Vfj) has a maximum value of 5.4 x 10 -14 at 
the mass upper limit of m(v^) = 0.19 MeV/c 2 Q. It 
is noteworthy that a BR of w 1CP 14 for 7r° — > vvj within 
the Standard Model is allowed. Therefore, an observed 
BR B(tt° — > UpV/j) » 5 x 10~ 14 would imply physics 
beyond the Standard Model. 

To date limits have been set on ir° — > vv derived from 
pion production in beam stops. An experimental upper 
limit, r(7r° -> Vyfip) /r(7r° -> all) < 3.1 x 10" 6 at 90% 
confidence level (CL), was set by Hoffman by using 
the data from several beam-dump experiments. Similar 
limits were obtained by Dorenbosch et al. ■ An inclu- 
sive search for ir° — ► vv using K + — > 7r + 7r° has set an 
upper limit of 8.3 x 10~ 7 (90% CL) 0. 

The Liquid Scintillator Neutrino Detector (LSND) ex- 
periment was performed at the Los Alamos Neutron Scat- 
tering Center (LANSCE). An 800-MeV beam of protons 
incident on a predominantly water-copper target pro- 
vided the source for the ir°. This energy is well below 
kaon production threshold, and almost all pions are pro- 



duced from p-nucleus interactions. The detector, located 
30 m from the beam stop, contained 167 tons of dilute liq- 
uid scintillator which served as the active target. The liq- 
uid scintillator was viewed by 1220, 8"-diameter Hama- 
matsu photomultiplier tubes (PMTs) mounted inside the 
tank. An active shield with 292, 5" -diameter PMTs ve- 
toed cosmic rays Q . The read-out and the data acquisi- 
tion system are described in more detail elsewhere p. 

The 7T° production rate was calculated using the LSND 
beam Monte Carlo (MC) Q. As the 7r° spectrum is ex- 
pected to be very similar to that of the ir + , in these 
calculations the tt + — ► e + v e decay mode was used, where 
a lifetime of 8.4 x 10~ 16 s (7r° lifetime) was assigned to the 
7r + . Assuming that all 7r°'s decay to two neutrinos, the 
neutrino flux is calculated to be (6.5 ±1.6) x 10 14 ^ M /cm 2 
at the center of the detector above the muon production 
threshold of 123 MeV for the entire running time of the 
LSND detector. The error estimate is based on the val- 
ues quoted by Burman and Plischke 0. Fig. ^shows the 
neutrino spectrum resulting from the possible n° — > v^v^ 
decay at the center of the LSND detector. Note that the 
peak between 70 and 80 MeV is due to multiple pion 
production 0- For comparison, Fig. Q also shows the 
v^ spectrum from the usual n + — » decay-in-flight 
(DIF) multiplied by 10. 

The results presented here were obtained during the 
1993-1998 LSND operation periods. A total of 28,896 
Coulombs of protons were incident on the beam stop. 
The dominant interactions inside the LSND detector 
were ^+ 12 C-> ^ +p + X, Z/^+ 12 C -> fi + +n + X and 
Vfi + p — > + n. The flux- averaged cross sections for 
these reactions are estimated to be (3.0±0.6) x 10 _38 cm 2 , 
(0.60 ± 0.12) xl0~ 38 cm 2 , and (0.20 ± 0.01) x 10~ 38 cm 2 , 
respectively, from a Fermi-Gas model (FGM) of Gaisser 
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TABLE I: Selection criteria and corresponding efficiency for 
the primary Michel electron. 



Selection Criteria 




Efficiency 


Veto Live Time 




0.76 ±0.02 


DAQ Live Time 




0.96 ±0.02 


Analysis Efficiency 


Electron Energy 


20 < E < 55 MeV 0.9 


8 ±0.02 


Fiducial Volume 


D > 35 cm 


0.8 


8 ±0.02 


Particle ID 


-1.5 < Xi < 0.5 


0.8 


4±0.01 


Shield and Crack Hits 


< 4 


0.9 


8±0.01 


Future Time Gate 


f uture 9 fiS 


0.99 ±0.01 


In-time Veto Gate 


Attto > 30 ns 


0.97 ±0.01 


Total in-time Cuts 




0.50 ±0.03 



TABLE II: Selection criteria and corresponding efficiency for 
the parent muon. 



FIG. 1: The neutrino energy spectrum resulting from the 
possible 7T° — » VpVp decay at the beam stop and as seen at 
the center of the LSND detector. The dashed curve shows the 
energy spectrum of neutrinos from ir + DIF. 



and O'Connell 9]. The estimated errors are from the 
uncertainties for the FGM calculations quoted by Vogel 
(Ipf . The n in each of the above reactions usually decays 
and produces a Michel electron. The 7r° — > v^v^ search 
focused on identifying high-energy, muon-like beam ex- 
cess events in the energy range between 160 MeV and 
600 MeV electron equivalent (MeVee). The lower cut of 
160 MeVee is obtained because the end-point energy of 
muon neutrinos is about 280 MeV and the muon pro- 
duction threshold is 123 MeV. This cut, therefore, in- 
sures that beam-related background /i events are negli- 
gible 0| • The upper energy cut was chosen because the 
end-point energy of the neutrinos is about 730 MeV and 
the muon production threshold is 123 MeV. 

The majority of the triggers (99%) were cosmic-ray in- 
duced and were classified as through-going muons, elec- 
trons from stopping muon decays, recoil protons from 
neutron collisions, and 12 i? g s . decay resulting from /i~ 
capture on 12 C. In order to minimize the background 
induced by cosmic rays, a Michel electron was selected 
for the current event that satisfied an energy cut of 
20 < E < 55 MeV. The reconstructed electron vertex 
was required to be inside a fiducial volume 35 cm from 
the face of the PMTs. The particle identification param- 
eter, Xii from reference |l2l | was used to select electrons 
and was required to satisfy —1.5 < xi < 0-5: where the 
allowed range is chosen by maximizing the selection effi- 
ciency divided by the square root of the beam-off back- 



Selection Criteria 



Efficiency 



Muon decay time 
Uncaptured pT 
Past Energy 
Spatial Correlation 
Fiducial Volume 



0.7 < t m < 9 lis 0.74±0.01 

0.97±0.01 

160 MeV < E < 600 MeV 0.64±0.01 
< 100 cm 0.99±0.01 
> 35 cm 0.75±0.01 



Shield and Crack Hits < 2 



0.87±0.01 



Total Past Cuts 



0.30±0.01 



ground. A veto cut of less than 4 hits was applied to re- 
duce cosmic rays. Events with future activity within 9 /is 
4// lifetimes) or with a bottom veto counter hit were 
rejected in order to further eliminate cosmic-ray muon 
events. Furthermore, no veto hit was allowed within 30 
ns of the trigger time. These cuts and their associated 
efficiencies for the selection of primary electrons are sum- 
marized in Table |U A second set of cuts was applied to 
isolate parent muons. These cuts are similar to those of 
reference [11]. Muons were required to have a decay time 
of less than 9/is. An energy cut between 160 MeV and 
600 MeV was chosen to eliminate possible muons from 
beam- related tt + DIF. A spatial [i — e correlation dis- 
tance cut of less than 100 cm was applied. The muon 
vertex was required to be within a fiducial volume of 35 
cm from the face of the PMT. A veto cut of less than 
2 hits was applied to further reduce cosmic rays. These 
cuts and their corresponding efficiencies are summarized 
in Table HD 

The overall efficiency of these cuts was 0.15 ± 0.01, 
which resulted in a total of 38 beam-on events and 473 
beam-off events. After applying the duty factor ratio 
of 0.062 ± 0.005, the normalized number of beam-off 
events was 29.3 ± 1.3 ± 2.4, resulting in a beam-excess of 
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FIG. 2: The visible energy distribution for beam-on, beam- 
off, and beam-excess events. The solid line represents the MC 
calculation. 



FIG. 3: Distance between muon and electron for beam-on, 
beam-off and beam-excess events. The solid line represents 
the MC calculation. 



8.7±6.3±2.4 events. Fig. [21 shows the visible energy dis- 
tribution for beam-on, beam-off, and beam-excess events. 
Also shown in the beam-excess figure is the area normal- 
ized GEANT-3.21 MC calculation. The calculation uses 
the neutrino spectrum of Fig. 1 and generates muons 
and knock-out protons based on the Fermi-gas model of 
reference [9]. 

Fig. |3|shows the spatial correlation between the muons 
and the Michel electrons for the beam-off, beam-on, and 
beam-excess events. The MC calculation is shown as 
the solid line on the beam-excess histogram. Fig. 0] 
shows the time correlation between the muons and the 
Michel electrons for the beam-on, the beam-off and the 
beam-excess events. The usual decay-at-rest muon life- 
time curves, e '/ 212Als are superimposed for comparison. 
Note the choice of 2.12 /is for muon lifetime is because 
the detected muons are a mixture of pT and fi + . 



The upper limit to the BR for the 7r° 



decay 



is calculated using the values for the neutrino reaction 
cross section (3.8 ± 0.8) x 10 _38 cm 2 ), the neutrino flux 
((6.5 ± 1.6) x 10 14 i//cm 2 ), the total number of target 
atoms ((3.7±0.1)xl0 30 ), and the overall efficiency (0.15± 
0.01). If the small excess is due to tt° — > v^v^ decay, 
then the corresponding BR is T (tt° ^v^p^) /T (tt° ^all) = 
[6A±4.6(stat)±3.3(syst)] x 10~ 7 , which corresponds to 
an upper limit of Ti^-^iy^D^/T^^all) < 1.6 x 10~ 6 
at 90% CL. 

In summary, a beam-excess of 8.7 ± 6.3 (stat) ±2.4 
(syst) events is observed above an energy of 160 MeV, 



where no beam related background is expected. Within 
the framework of ir° — ► vv, a direct upper limit for this 
BR of 1.6 x 10~ 6 at 90% CL is obtained. 
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FIG. 4: Muon lifetime for beam-on, beam-off and beam-excess 
events. The superimposed solid curves represent the expected 
e -t/2.l2fis ex p 0nen ti a l decay of muons at rest. 



